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Abstract. We show how the investigation of the redshifted 21-cm radiation can 
give insight into the development of structures in the early universe (at redshifts 
z > 5). In particular we investigate: the epoch of the first light; the fluctuations 
in the redshifted 21— cm emission induced by the density inhomogeneities in CDM 
dominated universes; the emission and absorption shells that are generated around 
the first bright quasars. Such features can be observed with the next generation radio 
facilities. 



1. The general framework 

The diffuse Intergalactic Medium (IGM) at very high redshift (between recombination 
and full reionization at z > 5) can be observed in the redshifted 21-cm radiation 
against the cosmic background. The signal can be detected in emission or in absorption 
depending on whether the spin temperature Ts is larger or smaller than the cosmic 
background temperature Tcmb = 2.73 (1 + z). 

This can happen if Ts is coupled via collisions to the kinetic temperature Tk of the 
IGM. However, the density contrast on Mpc scales at very early epochs is so low that 
the collision coupling is inefficient, and therefore the IGM is expected to be invisible 
against the CMB (Madau, Meiksin & Rees, 1998 lljl, hereafter MMR). On the other 
hand, large, massive regions at high density contrast are extremely rare in most of the 
hierarchical CDM universes at such high redshifts. 

There is a however another mechanism that makes the diffuse hydrogen visible in 
the redshifted 21cm line: the Wouthuysen-Field effect. In this process, a Lya photon 
field mixes the hyperfine levels of neutral hydrogen in its ground state via intermediate 
transition to the 2p state. A detailed picture of the Wouthuysen-Field effect can be 
found in Meiksin (1999, |§) and Tozzi, Madau, Meiksin & Rees (1999, [§, hereafter 
TMMR). The process effectively couples Ts to the color temperature T^ of a given 
Lya radiation field (Field 1958, B). The color temperature is easily driven toward the 
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kinetic temperature of the diffuse IGM due to the large cross section for resonant 
scattering (Field 1959, Q). In this case the spin temperature is: 

i + y. ' 

where j/q — 3.6 lO^^Pa/TK , and Pa is the total rate at which Lya photons are scattered 
by an hydrogen atom. 

However, the same Lya photon field also re-heats the diffuse gas, driving Tk 
toward larger values. The thermal history of the diffuse IGM then results from the 
competition between adiabatic cooling due to the cosmic expansion and re-heating due 
to the photon field. In the absence of a contribution from a strong X-ray background, 
the thermal history of the IGM can be written simply as: 

dTK _2^i E dt Tk , . 

dz 3 kgdz {1 + z)' ^ ' 

where E is the heating rate due to recoil of scattered Lya photons. Here /i = 16/13 
is the mean molecular weight for a neutral gas with a fractional abundance by mass 
of hydrogen equal to 0.75. 

Prior to the generation of the photon field, the IGM is neutral and cold, at a 
temperature Tk — 2.6 10~^(1 + z)^ (Couchman 1985, [^) given only by the adiabatic 
cooling after recombination. At the onset of the re-heating sources, there will be 
coupling between the kinetic and the spin temperature. An observation at the 
frequency 1420/(1 + z) MHz will detect absorption or emission against the CMB, 
with a variation in brightness temperature with respect to the CMB value: 

AT, . ,2..„,K) .-S(fg) ^ ^ . (3) 

where flu = ^ — — is the curvature contribution to the present density 
parameter, fl\ is the cosmological constant, fib is the baryon density, and rj = 
{TcMB - Ts)/Ts. 

Observations of such variation in the brightness temperature can be used to 
investigate the thermal history of the IGM, and the underlying birth and evolution of 
the radiation sources. All the following results are presented and discussed in TMMR 
(§)■ 



2. The epoch of the First Light 

We first investigate a simple situation in which, at a given redshift Zth, the Lya photon 
field reaches a thermalization rate Pth ~ 7.6 x 10"^'^ s~^ (1 + z); for such a value Tg 
is driven effectively toward Tk (see MMR, [Q). If the IGM is heated only by the 
same Lya photons, there will be a transient epoch where Ts < Tk, (i.e., 77 < 0) 
and an absorption feature appears at the corresponding redshifted frequency. This 
effect necessarily has a limited extension in time and thus in frequency space, since 
Ts becomes larger than Tqmb on a relatively short timescale. However, the signal 
is easily detectable with a resolution of few MHz, and, most of all, has a large effect 
since |?7| >> 1 when Ts « Tcmb- Such a strong feature marks the transition from 
a cold and dark universe, to a universe populated with radiation sources. 

If we assume that the Lya field reaches the thermalization value when 04;^ = 9 on 
a timescale r ~ 10 Myrs, the IGM will be visible in absorption for ~ 10 30 Myrs, 
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Figure 1. Top panel: ATj, for a resolution of 1' and 1 MHz in frequency, both for 
the fluctuations (dashed) and for the continuous distribution of the IGM (continuous 
line), assuming Pa to be equal to the thermalization rate at z^f^ 2; 9 in a critical 
universe. Bottom panel: the corresponding Ts (continuous) and Tk (dashed) are 
shown together with Tcmb (dotted line). 



corresponding to AT;, ~ 40 mK over a range of ~ 5 MHz. In the top panel of figure 
|l| such a signature is shown as a function of the observed frequency. In the bottom 
panel the corresponding thermal evolution for the IGM is shown. 

Such results are weakly dependent on the epoch Zth and on the adopted cosmology. 
However, the amplitude of the detected signal will be strongly dependent on the 
timescale r on which the Lya field reaches the thermalization value. In figure |^ the 
maximum of the absorption is plotted for different timescales r in three representative 
cosmologies. The signal is always larger than 10 mK; note however that for t > 30 
Myr, the absorption is spread out over a large interval in frequency, especially at 
v < 100 MHz where the sensitivity of radio telescopes becomes lower (see TMMR). 



3. The density field 

After re-heating and before reionization, Ts >> Tcmb holds, and the IGM is 
detectable only in emission. However, "q < \ always, and the effect due to the 
continuum distribution of a diffuse IGM is not as strong as in the absorption case. 
Such a small positive offset with respect to the CMB background can be difficult to 
detect. On the other hand, fluctuations in the redshifted 21-cm emission, which reflect 
fluctuations in the density of the IGM, are at least two orders of magnitude larger 
than the intrinsic CMB fluctuations on scales of ~ 1 arcmin. 

These fluctuations correspond to scales of a few comoving Mpc, and are in the 
linear regime at 2: > 5. In this case the fluctuations induced in the brightness 
temperature will be directly proportional to Ap/p, allowing a straightforward 
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Figure 2. The maximum of the absorption is shown as a function of the timescale 
on which the Lya photon field reaches (exponentially) the thermalization value at 
the same z^j, = 9. The different curves corresponds to Qq = 1, h = 0.5 (continuous 
line), Qo = 0.4, h = 0.65 (dashed line), and Uq = 0.3, A = 0.7, h = 0.7 (dotted line). 
The vertical lines mark the age of the universe at z = 9 in the three cases. 



reconstruction of the perturbation field at that epoch. In figure ^ and ^ we show 
results for two cosmologies, a tilted CDM universe with critical density (tCDM), and 
an open = 0.4 universe (OCDM). In both cases the fluctuations are normalized 
to reproduce the local abundance of clusters of galaxies. In OCDM the fluctuations 
are much larger (a factor of 3) since the evolution of the perturbation is strongly 
suppressed in an open universe with respect to the critical case, and for a given 
local normalization, the amplitude of the perturbations at high z is correspondingly 
larger. In both figures the density field has been evolved with a collisionless N-body 
simulation of 64^ particles using the Hydra code (Couchman, Thomas, & Pearce 1995, 
0). The box size is 20h~^ comoving Mpc, corresponding to 17 (11) arcmin in tCDM 
(OCDM). The baryons are assumed to trace the dark matter distribution without any 
biasing. Since the level of fiuctuations ranges from a few to ~ 10 /iJy per beam (with a 
resolution of 2 arcmin), it seems possible that observations with the Square Kilometer 
Array (Braun 1998, see also http: / / www.nfra.nl.skai ) may be used to reconstruct 
the matter density field at redshifts between the epoch probed by galaxy surveys and 
recombination, on scales as small as 0.5 — 2 h^^ comoving Mpc, i.e. masses in the 
range between 10^^ and 10^^ H'^Mq. 



4. The first quasars 

If re-heating is provided by a single quasar (without any other source of radiation) , 
21-cm emission on Mpc scales will be produced in the quasar neighborhood (outside 
the HII bubble) as the medium surrounding it is heated to Ts — Tk > Icmb by 
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Figure 3. Radio map of redshifted 21-cm emission against the CMB in a tCDM 
cosmology at ^ = 8.5. The linear size of the box is 20h~^ (comoving) Mpc. The point 
spread function of the synthesized beam is assumed to be a spherical top— hat with 
a width of 2 arcmin. The frequency window is 1 MHz around a central frequency 
of 150 MHz. The color intensity goes from 1 to 6 fi3y per beam. For clarity, the 
contour levels outline regions with signal greater than 4 ^Jy per beam. 
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Figure 4. Same as figure m for OCDM. 
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Quasar turning onatz = 9 

Total ionization rate of l(f photons/sec, spectral index of 1.8 
Critical Universe, h=0.5, n^ =0.8 
Box size 100 h' Mpc comoving 



z = 8.842 z = 8.267 z = 7.654 

lOMyr 50 Myr 100 Myr 




Figure 5. 21-cm emission against the CMB from the region surrounding a quasar 
source (lower left corner, in the center of an HII zone), revealed once the IGM is 
heated above the CMB by soft X-rays from the quasar. The angular resolution is 2 
arcmin, the frequency depth is IMHz, and the color levels range from to 10 mK 
with respect to the CMB level (dark). The temperature of the IGM is assumed to 
be Tff = Tqmb '^^ large distances from the quasar. 

soft X-rays from the quasar itself. The size and intensity of the detectable 21-cni 
region will depend on the quasar luminosity and age. In particular the intensity in 
the emission weakens with radius and with the age of the quasar. 

We calculated the kinetic temperature around a typical quasar, along with the 
neutral IGM fraction and Lya flux. The resulting radial temperature profiles were 
then superimposed on the surrounding density fluctuations as computed using Hydra. 
In figure ^ a sequence of snapshots after 10, 50 and 100 Myr after the birth of a quasar 
at z = 8.5 are shown in a box of lOOft."^ Mpc (comoving). The visual effect is due to 
the convolution of the spin temperature profile with the (linearly) perturbed density 
field around the quasar. The temperature of the IGM at great distances from the 
quasar is assumed to be Tk — Tcmb, and the signal goes to zero. In the figure the 
color ranges from to 10 mK with respect to the CMB level, which is black. 

Another situation occurs when the temperature of the IGM at large distances 
from the quasar is lower than Tcmb, e.g., Tk = 2.6 10~^(1 -I- z)^ K. In this case the 
emission region is followed by an absorption ring, since the Lya photons reach regions 
where Tk < Tcmb- The radio map resulting from a quasar 'sphere of influence' 10 
Myr after it turns on at z = 8.5 (tCDM) is shown in figure |[ The signal ranges 
from about —3 /iJy to 3/LtJy per beam (with a 2 arcmin resolution). The absorption 
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Figure 6. 21-cm emission and absorption against the CMB from the region 
surrounding a quasar source (lower right corner, in the center of an HII zone), 
revealed once the IGM is heated above the CMB by soft X-rays from the quasar. 
The linear size of the box is of 20h~^ comoving Mpc (tCDM), the angular resolution 
is 2 arcmin, the frequency depth 1 MHz. The color levels range from —3 fi3y to 
3 /ijy per beam. The contours mark the absorption ring. The quasar turns on at 
z = 8.5 with a ionizing photon luminosity of 10^^ photons s^^, and is observed after 
10 Myr. The temperature of the IGM beyond the light radius is assumed to be 
Tk ~ 2.610-2(1 + 2)2 ~ 2.6 K < TcMB- 
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Figure 7. Same as figure but the spectrum of the quasar has an exponential 
cutoff at energies larger than the Lyman limit. The absorption ring is much more 
evident and deep. 
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region is limited to a very sharp edge. However, in figure |j we show a quasar with 
the same Lya luminosity but with an intrinsic exponentially absorbed spectrum at 
energies larger than the Lyman limit. Consequently the HII region is reduced, and 
the X-ray warming front is well behind the light radius. This occurrence leads to a 
larger absorption ring where the signal reaches ~ — 20 /xJy in a 2 arcmin beam. 

Imaging the gas surrounding a quasar in 21-cm emission could provide a direct 
means of measuring intrinsic properties of the source, like the emitted spectrum and 
the opening angle of quasar emission. All these features are within reach of the new 
generation radio telescopes like SKA. 



5. Conclusions 



The Wouthuysen-Field effect allows one to peer into the Dark Age. The observation 
of the neutral IGM in the redshifted 21-cm can give insight into the thermal evolution 
of the diffuse hydrogen and thus into the formation and evolution of radiation sources, 
at epochs when the age of the universe is only ~ 0.3 Gyrs. In particular, the epoch of 
the First Light can be seen as a deep (~ 40 mK) absorption feature a few MHz wide 
against the CMB, at the corresponding redshifted 21cm line. Moreover, the density 
perturbation field at a redshift z « 5 20 can be reconstructed looking for mK 
fluctuations at 1-5 arcmin resolution in the radio sky, providing a determination of its 
amplitude between the epoch probed by galaxy surveys and recombination. Finally, 
the first ionizing sources, like luminous quasars, can be seen by identifying peculiar, 
ring-shaped signals whose morphology depends on the source's age, luminosity and 
geometry. 
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